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Abstract 

We discuss the phenomenon of spontaneous symmetry breaking by means of 
a class of non-perturbative renormalization group flow equations which employ a 
regulating smearing function in the proper-time integration. We show, both an- 
alytically and numerically, that the convexity property of the renormalized local 
potential is obtained by means of the integration of arbitrarily low momenta in 
the flow equation. Hybrid Monte Carlo simulations are performed to compare the 
lattice Effective Potential with the numerical solution of the renormalization group 
flow equation. We find very good agreement both in the strong and in the weak 
coupling regime. 



1 Introduction 



The mechanism of spontaneous symmetry breaking is one of the most important nonper- 
turbative phenomena in quantum field theory and statistical mechanics. Unfortunately, 
due to its strong nonperturbative character, many of its aspects are not yet well under- 
stood mainly because of the lack of reliable computational tools even in simple physical 
situations. For example, in the liquid-vapor phase transition, the system undergoes a 
first-order phase transition below the critical temperature, but the presence of long-range 
fluctuations renders the mean-field description inadequate. The well-known Maxwell con- 
struction if] reproduces the fiat isotherms in the two-phase region but does not provide 
the correct critical exponents near the critical region. 

In quantum field theory the situation is even more dramatic because the perturbative 
calculation of the Effective Potential (EP) generates an imaginary part whose physical 
meaning is not entirely clear. In fact, the EP should be a convex function of the fields 
by construction, since it is the Legendre transform of the generating functional of the 
connected Green's functions j2j. 

An important result was obtained in ^ within the framework of the Effective Average 
Action, for a continuous 0(A^)-symmetric scalar theory: it was shown that in order to 
correctly reproduce the convexity property of the free energy a non-trivial saddle point 
(spin wave solution) must be chosen in the loop-expansion P^. This result has opened the 
possibility of using the Wilsonian RG approach in discussing non-perturbative issues of 
this kind. 

Quite generally the possibility of applying renormalization group (RG) methods to 
the description of first order phase transition has frequently been matter of discussion in 
the recent years. Although the question is still not completely settled, several works have 
shown that the Wilsonian RG approach is a promising tool for such kind of investigations. 
In particular it has been shown [S^ that the Wegner-Houghton (WH) RG equation (a sharp 
cut-off momentum-space wilsonian-type of RG flow equation) already in the simple local 
potential approximation (LPA) reproduces the convexity property of the free energy as 
a result of the integration of arbitrarily low momenta without resorting to any ad hoc 
Maxwell construction, or saddle point approximation. Moreover, the RG transformation 
is always well defined below the critical line. 



1 



The essential ingredient of their analysis was the use of a robust and very accurate 
numerical method for handling the RG equation without resorting to polynomial trun- 
cations. Similar results have also been obtained with a smooth cutoff RG flow equation, 
in the framework of the exact RG equation for the effective average action 6J. In [7 an 
alternative method which combines the RG flow equation approach with the non-trivial 
saddle point expansion method has been presented. 

The wilsonian RG transformation [9^ , as opposed to the more conventional RG trans- 
formation based on the rescaling properties of the Green's functions, preserves all the 
infinite number of interactions generated in the low-energy action by "averaging-out" the 
small-scale physics defined at some high-energy scale A, and for this reason is a power- 
ful tool in order to investigate situations where infinitely many length scales are coupled 
together. Quite generally, the Wilsonian action can be defined as 

e-^"(*) = J Vip5{^ - C{ip)) e-^(^) (1.1) 

where $ is a low frequency field which is in turn a function of the fundamental field 
If via $ = C(0) and C is some averaging operator. If the effective theory is defined on a 
lattice Aq, then 

For computational purposes, a continuous wilsonian RG transformation is defined in 
the momentum-space representation where the definition p.lj) becomes 

e-5.(*) ^ J p^e-^K+*) (1.2) 

( being the high-frequency part of the fundamental field with momenta p greater than 
some infrared cutoff p > k, and $ the low-frequency one, with momenta p < k. If we 
choose 

Civ) = ^ f ^ (1-3) 

where is a characteristic volume over which the field is averaged, in the limit Q oo 
(equivalently, vanishing infrared cutoff k in Eq. ()1.2|1 ) the blocked action gives precisely the 
EP, namely the non-derivative part of the Effective Action r[0], defined as the Legendre 
transform of Vr[J], the generator of the connected Green's functions in the presence of 
an external current J jSI- 



2 



In practice, the wilsonian RG transformation obtained from explicit calculation of 
(ll.2|l is strongly non-linear, and only by retaining infinitely many interactions generated 
by integrating out the "fast" variable, the RG transformation can be well defined below 
the critical line. Therefore, a numerical treatment of the full non-linear partial differential 
equation is necessary even in the simplest LPA approximation of the Wilsonian Action. 

In recent years, in an attempt of calculating ()1.2|1 non-perturbatively, a new type of 
RG flow equations has been proposed [?] and further developed jUIl HH [121 • based on 
Schwinger's proper time regulator, and it is therefore well suited for a direct application 
to gauge theories. 

The proper time renormalization group equation (PTRG) does not belong to the class 
of exact-RG flows which can be formally derived from the Green's functions generator 
without resorting to any truncation or approximation f!M [T^ [W\ . Nevertheless the flow 
equation has the property of preserving the symmetries of the theory and, in addition, 
it has a quite simple and manageable structure. Moreover the PTRG flow, although 
it can not fully reproduce perturbative expansions beyond the one loop order [T^, still 
does provide excellent results when used to evaluate the critical properties such as the 
critical exponents of the three dimensional scalar theories at the non-gaussian flxed point 
^HIIEI- In particular the determination is optimized when one takes the "sharp" limit 
of the cutoff function on the proper time and it turns out to be much more accurate 
than the one corresponding to a smoother regulator. The PTRG also provides excellent 
determinations of the energy levels of the quantum double well j2^. In fact, the use of 
the special type of PT regulator described in this paper and already discussed in [iHl EHI 
realizes an localized integration over the "fast" variables at a given infrared running 
momentum k and, for this reason, it shows an impressive numerical convergence stability 
in the calculations of the critical exponents. 

In this paper we shall instead be mainly concerned with the study of a a strong non- 
perturbative phenomenon like the SSB, for which the validity of PTRG equations has 
not been tested yet. In particular, the relevant question we would like to address in this 
paper is to analytically and numerically study the approach to the convexity for the EP 
obtained by means of a class of PTRG flow equations in the limit of vanishing infrared 
cutoff. A smeared type of proper-time regulator will be used which selects a localized 
integration in the "fast" variables in the path-integral. The resulting flow equation can 
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be thought of as an interpolation of a "smooth" modification of the LPA approximation of 
the Wegner- Houghton (WH) [12] "exact" flow equation, and the exponentially ultraviolet 
(UV) converging flow discussed in ^[Ej. It will be shown that the PTRGs can correctly 
describe the approach to convexity as A; ^ 0, and that give results which are completely 
consistent with the analysis performed in for the WH equations, and in j3] for the 
"exact" evolution equation for the effective average action. The additional advantage is 
that below the upper critical dimension the PTRG flow equation predicts the expected 
discontinuity of the inverse compressibility, whereas the "exact" WH equation fails to do 
so. 

In order to strengthen our analysis, we compare the result of our numerical investiga- 
tion with the standard EP computed on a lattice by means of a Hybrid Monte Carlo ap- 
proach. Our main conclusion is that lattice EP and PTRG flow equation agree extremely 
well even for relatively small lattices, both in the perturbative and non-perturbative 
regime. 



2 Flow equations from proper-time regulators 

Let us now review the basic assumptions in the derivation of the PTRG flow equation. 
We first calculate the wilsonian action in ()1.2|) in the one-loop approximation as 



6^S 



Sl-''""m = -^ ^/,(.,A)Tr(e ^ - e ^ (2.1) 

where, following [211 122 I2H| , we have introduced an heat-kernel smooth regulator fk 
in the Schwinger proper time representation for the one-loop expression and we have 
subtracted the UV divergent contribution of a vanishing background field. The following 
prescription apply to the choice of fk'- 

• /fc(s = 0, A) = in order to regularize the UV divergences as s — > 

• fk=o{s ^ oo. A) = 1 in order to leave the IR physics unchanged by the regulator 

• fk=A{s, A) = in order to recover the original bare theory at the cutoff scale k = A. 

The flow equation for the Wilsonian action 5*^ results from considering Eq. ()2.1|) in the 
infinitesimal momentum shell k,k + Sk and then RG-improving the resulting expression, 
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which amounts to substitute the bare action S\ with the running one S\ Sk in the 
RHS of Eq.dSH), 

dk 2 Jo s dk \ J ^ ^ 

where primes indicate functional derivatives with respect to (p. At this point, one 
should stress that although the fast momentum integration in ()2.H1 has been obtained by 
assuming a trivial saddle point in the path-integral, the RG-improved flow equation ()2.2j) 
is instead a functional differential flow equation which is valid for a general average action 

If we then project Eq. fl2.2j] onto its non- derivative part, the anomalous dimension rj is 
set to zero, and the resulting flow equation for the local potential reads 

dUkm K, p ds ,dMs,A) 



k- 



dk 







where Kd = 2/(47r)'^/^r(c?/2). For actual calculations we shall employ the following class 
of proper-time smearing functions 

;.(,A)^£(!!ii^!y(!^ ,2,4) 

where r(a, z) = dt t"~^e~* is the incomplete Gamma function and n a positive real 
number. Note that lim^^oo = ('^s. A) = e{s - k'^) - 9{s - A^). 

Since kdT{n, sk'^)/dk = — 2(sA;^)"e~'^'^'^ the s-integration in ()2.3|1 can be explicitly 
performed, having at last 

where = n^^'^T{n — d/2) /T(n), and in particular = 1. In order to simplify the 
notation we have not explicitly written in the RHS of the equations the field independent 
contributions which correspond to a trivial shift of the vacuum energy at /c = 0. Eq. ()2.5p 
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is obtained from the n = d/2 limit by performing the trivial rescaling k ^Jd/2k of the 
infrared cutoff. It coincides with the LPA approximation of the "exact" WH sharp cut-off 
equation Eq. ()2.(i|l converges to Eq. ()2.7|1 as n ^ oo. In other words, the n- dependence 
in the cutoff interpolates between a WH type of equation and the "exponential" RG 
equation ()2.7|) . 

3 Analytical and numerical results 

We shall now analytically and numerically study the RG flow below the critical line. 
Eq. ()2.2j) is in principle an evolution equation for a generic Wilsonian action Sk{^) and 
already in the simplest LPA approximation described in ()2.5ll2.7p it directly exhibits the 
approach to convexity as A; ^ 0. In other words, the non-perturbative features like spin 
waves, kinks or instantons are already included in Eq. ()2.2j) . This result should not come 
entirely as a surprise since we already stressed that ()2.5j) obtained for n = d/2, is the LPA 
approximation of the "exact" Wegner- Houghton equation studied in [5J. In this respect, 
Eq. ()2.6p and Eq. ()2.7|) can be thought of as smooth-cutoff modifications of the sharp cutoff 
WH equation, and it is then important to understand their IR behaviour below the critical 
temperature. In particular, the location of the polar singularity in Eq. ()2.(j|l in the ($, k"^) 
plane seems to be n-dependent: 

nk"" + = (3.1) 

so that in the limit n ^ oo, it becomes an essential singularity at = in Eq. ()2.7|) . The 
relevant question, is to understand how the convexity of the free energy is recovered in 
the — i> limit in these cases. 

Our discussion will be mainly concerned with n > d/2, because the n = d/2 case has 
extensively discussed by |5^^ and the n — oo case can be often extrapolated as a large n 
regime of ()2.6j) . (numerically one does not see any relevant difference between ()2.6j) and 
(jTTjl already for n ~ 20). 

In order to study the behavior of the solution of Eq. ()2.6|) near $ = it is convenient 
to introduce the rescaled potential and field 

Vk = M^/{47Tf^Uk = ^M^/(47r)^/2$. (3.2) 
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It is then convenient to take the second derivative of Eq. ()2.6|) with respect to and 
to define a new variable Wk{(p) through 

W,{<P) = 1 + -^) (3.3) 



which becomes large and positive in the broken phase near = for A; ^ 0. In terms of 
this new variable Eq. ()2.(jj) becomes 
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W'i = -2k-'^M^ - W;^') - ^^-^ W; "-^^ k^ (3.4) 

and prime means now derivative wrt (p. In Eq. ()3.4|) . when is large and positive, we 
can neglect terms which are suppressed as inverse power of Wk and obtain the solution 

W, = n k-'^'i^l - (3.5) 

with 101 < 00; which is well behaved for any value of A; > 0, being 0o an integration 
constant. Since our approach breaks down at = 0o, we can thus identify 0o as the value 
of the field at the coexistence. 

By inserting solution ()3.5|) back in ()3.3|) and solving for f/^' C^) have the well-known 
behavior 

f/fc(0) = -^0' + O(fcV) (3.6) 

which is consistent with the analysis of [3 H] . Incidentally, we note that ()3.6|) shows 
that the approach to the fiat bottom is slower for larger values of n. 

For the sake of completeness we notice that the n oo limit of ()3.4|1 reads 

<^2*-ln,y.-i^*^ (3.7) 

while if we instead define Wk{4>) = ln(l + Vl'/k"^), Eq. ()2.5p implies 

WH = 4fc-'^+2 _ 2e^^ k^. (3.8) 

Although no simple analytical solution is available in the broken phase near = for 
Eq. ()3.7j) we can still understand how the regulator affects the approach to the spinoidal 
line. 
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Figure 1: The blocked potential in ci = 3. The solid line is for n = 1.5 (WH), the dashed 
is for n = 2 and the dot-dashed is for n = 4 

In fact, when Wk is large and negative in Eq. ()3.8|) . the diffusive term dW/dk is expo- 
nentially suppressed in this case and we obtain Eq. ()3.(i|l with n = 1 (see also jS]). On the 
contrary when Wk is large and positive in Eq. ()3.4j) and Eq. ()3.7j) the diffusive term is only 
power-law suppressed. In particular the suppression is the slowest for — > cxd, like l/W^, 
while if n = d/2 + e (being e small and positive), the suppression is like W/. ^^'^^^^'^ and 
fl3.5|l is clearly reached faster as e gets smaller. In other words convexity is best achieved 
with a small value of e. 

An important issue related with the above discussion is to understand how the inner 
solution joins the outer region (p > (po. For Eq. ()2.5j) it is known that a fixed point solution 
for — > is present, but it predicts a continuous second derivative of the effective potential 
for d < A, which corresponds to a diverging compressibility at the transition (see [S] for 
an extended discussion on this point). The relevant question is whether the use of a 
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Figure 2: The blocked potential in d = A. The solid line is for n = 2 (WH), the dashed is 
for n = 3 and the dot-dashed is for n = 5. 

smooth cutoff regulator as Eq. ()2.6|) and Eq. ()2.7|) can cure this pathological behavior of 
the "exact" WH equation. 

In order to investigate these problems in detail one must handle the problem numer- 
ically, by extracting an accurate numerical solution of the flow equation. We thus have 
solved Eq. ()3.4|1 and Eq. ()3.7|) with the fully-implicit, predictor corrector finite-difference 
scheme described in |24 for which a rigorous result ensures convergence to the real solution 
for our numerical discretization grid. 

Let us then write the bare potential as Va(0) = to0^/2 + (y'o0'^/4!, and |ro| < 1 being 
the bare mass measured in cutoff units. Fig.(^ and Fig.Q show the blocked potential 
for (i = 3, and ci = 4, respectively, in the broken phase. As expected, if we push the 
integration closer to the A; — > limit all the curves approach a completely flat bottom as 
it is shown in Fig.(^ and Fig.(j2)) in = 3 and c? = 4, respectively, for various values of 
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Figure 3: The blocked potential in d = 3 for various values of the RG time t. 

n. In this case tq = —0.6, go = 1.0, and the final value of the RG "time" is t = —10. 
Different values of tq and Qq leads to qualitatively similar results in the sense that, as 
long as we are below the critical line, the flow always approaches the flat bottom convex 
potential. 

As we discussed above, we also notice from Fig.Q that the convergence to convexity 
is much faster for smaller values of n. A plot for different values of the RG time t is 
depicted in Fig. (jH)) where the expected behavior is discussed: as it is apparent from these 
plots the the potential for n = 2 is much closer to the fiat and convex solution already for 
t = 5, than the n = oo case where, although a fiat bottom is present near = 0, convexity 
is not achieved yet. It is numerically difficult to reach higher values of t while keeping the 
mesh spacing constant if n oo. We find that in order to reproduce the convexity at 
an acceptable level {\min{U) — U{0)\ ~ 10"^ for qq = 0(1)) the time step and the mesh 
spacing have to be both of the order of at least 10"'^, which is not very efficient from 
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Figure 4: V" as a function of (j) in the broken phase for d = 3 and with tq = —0.1 and 
go = 1.2. The squares are for n = 3 while the triangles are for the WH equation. The 
discontinuity is clearly visible in the first case, while it is not present for the WH equation. 

the numerical point of view. On the contrary, if we choose n = d/2 + e (e positive and 
small), then \min{U) — U{0)\ ~ 10^^ already with a mesh spacing one order of magnitude 
greater. Although we find numerical evidence that for n ^ oo the solution approaches 
a flat bottom near the origin, we cannot exclude that convexity is never reached in this 
case. 

An important result of our analysis is that the second derivative of the potential shows 
the expected discontinuity for 3 < c? < 4, for any n > d/2, as opposed to the "exact" WH 
equation, for which the spinoidal line merges with the coexistence line in d = 3. This is 
clearly shown in Fig.Q, where the discontinuity is visible in the numerical output because 
always only one grid point is present in the jump, and this feature does not depend on 
further refinements of the spatial grid. Similar behavior is also observed in d = 3.5, as it 
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Figure 5: V" as a function of in the broken phase for d = 3.5 and with tq = —0.1 
and go = 1.2. The squares are for n = 3 while the triangles are for the WH equation. 
The jump is clearly seen for the squares while there is a continuous transition for the WH 
equation 



is apparent from Fig.©. 



4 Monte Carlo simulations 

Monte Carlo simulations of the 0^ model have been performed on c? = 2, 3, 4 lattices in 
order to evaluate the Effective Potential and compare its determination with the numerical 
solution of the flow equations. In a lattice simulation, the inverse lattice spacing acts as 
an UV momentum cutoff, while the finite size introduces an IR cutoff. The model was 
discretised following the standard approach in and the field configurations were 
generated with the Hybrid Monte Carlo algorithm. Introducing an uniform external 
current J, a discrete version of the action can be written as follows: 
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Figure 6: The average field as evaluated by the HMC simulation, in the d = 3 case 
(L = 16), with parameters o^Tq = —0.2, ago = 0.24. The results for three different values 
of J are shown. 



1 ^ 

^ ^{'Pn+t, -~<PnY + V{(l)n) - J<Pr, 



where n is the site index, a is the lattice spacing and 



(4.1) 



V{<P)=a^^<P' + a^-'^-^<P^ (4.2) 

Introducing the normalization a'^/'^^^cj) (2k)^/^0 and a'^/^+^(2K)^/^ J J, the dis- 
cretised, dimensionless action is: 



where 



-2k J] <j)n+^.<Pn + 4>l + K4>1 - 1)' - A 



2 1 - 2A , 4_rf 6A 
aVo = 2d a" "go = ^ 



(4.3) 



(4.4) 



Periodic boundary conditions have been assumed. Both cold and hot starts have been 
tested, and once thermalization was achieved, the average of the field over the entire 
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Figure 7: The blocked potential from the flow equation in c? = 3 vs the MC simulation 
for a volume of 12^ in the weak coupling case (MC errors negligible on this scale). 

lattice was computed on each configuration: 

^-T^Y.^- (4-5) 

n 

where L is the number of lattice sites in each direction. The results were saved every 
50 configurations to decorrelate the results. In each case, 200 uncorrelated configurations 
were collected in order to get a statistically accurate determination of the average field. 
It should be noticed that the size of the statistical fluctuations decreased as the external 
current J was increased, as expected. In all cases, statistical errors have been estimated 
using the bootstrap technique (1000 samples) and found to be negligible. Simulations 
have been run with different lattice volumes, and flnite size effects have been found to 
be negligible. For each choice of the numerical parameters {k, A) , the simulation was run 
with several values of the external current J. This procedure was not very different from 
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Figure 8: The blocked potential from the flow equation in d = 3 vs the MC simulation 
for a volume of 12^ in the strong coupling case (MC error bars invisible on this scale). 

that sketched in 1^2] • The relation $(J) was then inverted to get and thus the first 

derivative of the Effective Potential by 

J = (4.6) 

We would like to compare U' as computed from ()4.6|) with the same quantity obtained 
by numerical integration of the flow equation. Although it would be possible to compare 
the constraint effective potential jHSI at a given scale k oc with the blocked potential 
as computed from the flow equation at a scale k, it is more interesting to focus on the 
A; ^ limit, where possible non- universal features due to the regulator may disappear. 
As we discussed in the introduction we prefer to use the external current method because 
it is simpler to implement and as accurate as the constraint effective potential method 
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We explore a set of parameters which is not close to the critical line, so that finite- 
size effects can be neglected, and we are able to compare the result of the fiow equation 
directly to the lattice determination of the effective potential. In particular we find that 
it is was not necessary to construct a lattice version of the RG fiow equation as discussed 
in [33]. We have integrated Eq. ()3.4j) rewriting all the relevant quantities in units of the 
UV cut-off and then we have followed the evolution down to t — > oo with the help of the 
numerical integrator. In order to show the predictive power of our approach we have not 
fine-tuned the bare parameters in the RG fiow equations to reproduce renormalized mass 
and coupling constant obtained in the lattice calculation as done in [HI] • Instead, we have 
decided to set the same bare mass and coupling constant in the lattice bare potential 
fl4.2|) and in the bare potential of the RG equation. Moreover we have rescaled back 
our potential and field according to ()3.2|) so that we get f/(.^Q($) out of the numerical 
computation. According to the analysis of the previous session we have considered n G 
{d/2,d\ where the numerical stability is best achieved. 

The results are shown in Fig.© and Fig.® for n = 2 and d = 3, where it is apparent 
that there is already a very good agreement with the MC data both in the weak and in 
the strong coupling regime. Better agreement could probably be achieved by including 
the wave-function renormalization function, but this is not our main concern in this 
investigation. 

5 Conclusions 

We have discussed the PTRG fiow equation below the critical line in a scalar theory. In 
particular we have shown that the convexity property of the free energy is recovered by 
integration of the LPA fiow equation in the k ^ oo limit. Within a class of n-dependent 
proper time regulator, the approach to the correct fiat bottom potential is faster when 
n = d/2 + e being e positive and small. The expected discontinuity of the U" at the 
transition is correctly reproduced for any value of e > as opposed to the "exact" WH 
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flow (e = 0), which does not show this feature. We have performed an extensive MC 
investigation of the EP in = 3 in order to discuss the numerical predictions and we 
found very good agreement both the strong coupling and weak coupling phase without 
resorting to a fine-tuning procedure between the bare parameters in the MC and in the 
RG flow equation. We anticipate that our result can be relevant in gauge theory where the 
presence of the PT regulator is an essential tool deriving a non-perturbative flow equation 
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